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Abstract. Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is simulated 
assuming ^So pairing of neutrons in the NS crust, and also ^So pairing of protons and weak ^P2 pairing of neutrons 
in the NS core, and using realistic density profiles of the superfluid critical temperatures T'c(p). The theoretical 
cooling models of isolated middle-aged NSs can be divided into three main types. (I) Low-mass, slowly cooling NSs 
where the direct Urea process of neutrino emission is either forbidden or almost fully suppressed by the proton 
superfluidity. (II) Medium-mass NSs which show moderate cooling via the direct Urea process suppressed by the 
proton superfluidity. (Ill) Massive NSs which show fast cooling via the direct Urea process weakly suppressed by 
superfluidity. Confronting the theory with observations we treat RX J0822-43, PSR 1055-52 and RX J1856-3754 
as slowly cooling NSs. To explain these sufficiently warm sources we need a density profile Tc{p) in the crust with 
a rather high and flat maximum and sharp wings. We treat IE 1207-52, RX J0002-I-62, PSR 0656+14, Vela, and 
Geminga as moderately cooling NSs. We can determine their masses for a given model of proton superfluidity, 
Tcp(p), and the equation of state in the NS core. No rapidly cooling NS has been observed so far. 
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1. Introduction 

Cooling of neutron stars (NSs) depends on the properties 
of matter of subnuclear and supranuclear density in the NS 
crusts and cores. These properties are still poorly known 
and cannot be predicted precisely by contemporary mi- 
croscopic theories. For instance, microscopic calculations 
of the equation of state (E OS) o f matter in the NS cores 
(e.g., Lattimer & Prakash ^001 ) or the superfluid prop- 
erties of NS cores and crusts (e.g., Lombardo & Schulze 
2001 ) show a large scatter of results depending on a model 
of strong interaction and a many-body theory employed. 
It is important that these properties can be studied by 
confronting the results of simulations of NS cooling with 
the observations of thermal emission from isolated middle- 
aged NSs. 

This paper is devoted to such studies. For simplicity, 
we use the NS models with the cores composed of neutrons 
(n) with an admixture of protons (p) and electrons. We 
will mainly focus on the superfluid properties of NS matter 
which are characterized by the density-dependent critical 
temperatures Tc(p) of nucleons. It is customary to consider 
superfluidities of three types: singlet-state (^Sq) superflu- 



idity (Tc = Tens) of neutrons in the inner NS crust and 
the outermost core; ^Sq proton superfluidity in the core 
{Tc = Tcp); and triplet-state ('^P2) neutron superfluidity in 
the core (Tc = Tent)- Superfluidity of nucleons suppresses 
neutrino processes involving nucleon s and affects nucleon 
heat capacities (e.g., Yakovlev et al. 1999). In addition, it 



initiates a specific mechanism of neutrino emission associ- 
ated with Cooper pairing of nucleons (Flowers et al. 1976 ). 
Our aim is to analyze which critical temperatures Tc(p) 
are consistent with observations and do not contradict the 
current microscopic calculations. 

We have considered this problem in two prior publica- 
tions. Kaminker et al. (2001, hereafter Paper I) analyzed 



the effects of proton superfluidity (basing on one partic- 
ular model of Tcp(p)) and '^P2 neutron superfluidity in 
the NS core. Yakovlev et al. (2001, hereafter Paper II) 
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included, additionally, the effects of So neutron superflu- 
idity in the crust. Calculations in Papers I and II were 
performed for one particular EOS in the NS core. In the 
present paper we extend the results of Papers I and II 
by considering three models of proton superfluidity and 
another EOS in the NS core. We combine the results of 
Papers I and II and give an overall analysis of the problem. 
We show that one can distinguish three distinctly differ- 
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ent types of isolated middle-aged NSs, which show slow, 
moderate, and fast cooling. Using this concept wc discuss a 
possible interpretation of observations of thermal emission 
from eight middle-aged isolated NSs. 

2. Cooling models 

We simulate NS cooling with our fully relativistic non- 
isothermal cooling code (Gnedin et al. 2001). The code 



Table 1. NS models employing EOSs A and B 



solves the radial heat diffusion equation in the NS interior 
(excluding the outer heat-blanketing layer placed at p < 
10^° g cm"'^). The heat is carried away by the neutrino 
emission from the entire stellar body and by the thermal 
photon emission from the surface. No additional reheating 
mechanisms are included. 

The code calculates theoretical cooling curves, i.e., the 
effective surface temperature as detected by a distant ob- 
server, T^, versus NS age t. The thermal history of an 
isolated NS consists of three stages. The first is the stage 
of thermal relaxation of the stellar interior (e.g., Lattimer 
et al. |1994 Gnedin et al. |200l|). It lasts for about 10-100 



yr. It is followed by the stage at which the NS interior is 
isothermal and the neutrino luminosity exceeds the sur- 
face photon luminosity (t ^ (3 — 10) x 10^ yr). The final 
stage is the photon cooling stage at which the photon lu- 
minosity dominates the neutrino one. 

In the NS crust we use the EOS of Negele & Vautherin 
1973D (atomic nuclei everywhere in the crust are assumed 



to be spherical). The core-crust interface is placed at the 
density 1.5 x 10^"* g cm~^. A standard procedure is used 
to match the core and crust EOSs near the core-crust in- 
terface. In the core, we use two p henomenological EOSs 
proposed by Prakash et al. ( 1988 ). We refer to them as 
EOS A and EOS B. 



EOS A is model I of Prakash et al. (|l988[) with the com- 
pression modulus of saturated nuclear matter K = 240 
MeV. It has been used in Papers I and II. EOS B corre- 
sponds to K = 180 MeV and to the simplified form of th e 
symmetry energy proposed by Page & Applegate ( 1992| ). 
EOS B has b een u sed in a number of pa pers (e.g.. Page 
& Applegate 1992| , Yakovlev et al. 1999 and references 
therein). 

The masses, central densities, and radii of two stellar 
configurations for each EOS are given in Table 1. The first 
configuration is the most massive stable NS. The values 
of Mmax indicate that EOS A is stiff while EOS B is mod- 
erate. The second configuration has a central density at 
which the direct Urea process switches on; for both EOSs 
it is allowed at M > Md (pc > Pd)- EOS B impHes a 
smaller symmetry energy at supranuclear densities and 
opens the direct Urea process at a higher density. 

Our cooling code includes all the important neutrino 
emission processes in the NS core (direct and modi- 
fied Urea, neutrino bremsstrahlung in nucleon-nucleon 
collisions, neutrino emission due to Cooper pairing of 
nucleons) and in the crust (plasmon decay, neutrino 
bremsstrahlung due to scattering of electrons off atomic 
nuclei, electron-positron annihilation into neutrino pairs. 



Model 


Main parameters 


EOS A 


EOS B 


Maximum 


Mmax/M© 


1.977 


1.73 


mass 


Pcmax/lO^'* g Cm"^ 


25.75 


32.5 


model 


R km 


10.754 


9.71 


Direct Urea 


Md/Mq 


1.358 


1.44 


threshold 


Pd/10" g cm-3 


7.851 


12.98 


model 


R km 


12.98 


11.54 



Table 2. Parameters of superfluid models in Eq. (1) 
To/lO^ K 



Pair- 
ing 



Mo- 
del 



'So 
'So 
'So 
'So 
'So 



Ip 

2p 

3p 

Ins 

2ns 

3ns 

Int 



20.29 

17 

14.5 

10.2 

7.9 

1800 

6.461 



fco 


fei 






fm'i 


fm-i 


fm-i 


fm-i 





1.117 


1.241 


0.1473 





1.117 


1.329 


0.1179 





1.117 


1.518 


0.1179 





0.6 


1.45 


0.1 





0.3 


1.45 


0.01 





21 


1.45 


0.4125 


1 


1.961 


2.755 


1.3 



neutrino emission due to Cooper pairing of free neutrons in 
the inner crust). The effects of nucleon superfluidity are in- 
corporated in the neutrino reaction rates and nucleon heat 
capacities as described in Yakovlev et al. ( 1999| , 2001). The 
effective masses of protons and neutrons in the core and 
free neutrons in the crust are taken to be 0.7 of the bare 
nucleon masses. The values of the thermal conductivity in 
the NS crust and core are the same as used by Gnedin et 
al. ( pOOll ). 

The relationship between the effective surface temper- 
ature and the temperature at the bottom of the outer 
heat-blanketing envelope is taken according to Potckhin 
et al. ( |l997l ) and Potekhin & Yakovlev ( |200l| ). This allows 
us to consider either the models of NSs with the surface 
layers made of iron, without magnetic field and with the 
dipole surface magnetic fields B ^ 10^^ G, or the non- 
magnetic NS models with the surface layers containing 
light elements. It is assumed that the surface magnetic 
field induces an anisotropic heat transport in the heat- 
blanketing envelope but does not violate the isotropic (ra- 
dial) heat diffusion in the deeper NS layers. It is also as- 
sumed that a NS may have a hydrogen atmosphere even if 
the heat-blanketing envelope is mostly made of iron. The 
majority of cooling curves will be calculated for the model 
of non-magnetized heat-blanketing envelope made of iron. 
Two exceptions are considered in Sect. 4.3. 

Although the microscopic calculations of superfluid 
critical temperatures Tc{p) give a large scatter of re- 
sults (Sect. 1), some common features are evident. For 
instance, Tc{p) increases with p at sufficiently low densi- 
ties (due to an increasing strength of the attractive part of 
nucleon-nucleon interaction), reaches maximum and then 
decreases (due to a short-range nucleon-nucleon repulsion) 
vanishing at a rather high density. For Tcnsif), the maxi- 
mum takes place at subnuclear densities, while the switch 
off occurs at p ^ po, where pa w 2.8 x 10^'* g cm~'^ is the 




Fig. 1. Density dependence of the critical temperatures 
for three models Ip, 2p, and 3p of the proton superfluidity 
(dots-and-dashes) in the core (with EOS A); three mod- 
els Ins, 2ns, and 3ns of ^So neutron superfluidity (solid, 
short-dashed, and long-dashcd lines); and one model Int 
of '^P2 neutron superfluidity (dots) used in cooling simu- 
lations. The parameters of the models are given in Table 
2. Vertical dotted lines indicate neutron drip point, core- 
crust interface, and the direct Urea threshold. 



Fig. 2. Superfluid gaps for ^Sq neutron pairing versus 
neutron Fermi wavenumber fcpn for models Ins, 2ns, and 
3ns (solid line Ins, short-dashed line 2ns, and long-dashed 
line 3ns; Fig. 1, Table 2). Sohd Hue 1 is the gap derived 
from the BCS theory with in- vacuum nn interaction (after 
Lombardo & Schulze 2001); dotted lines 2 and 3 are two 
different theoretical curves obtained by Wambach et al. 
( |1993| ) and Schulze et al. ( 1996 ) including medium polar- 
ization effects. 



saturated nuclear matter density. For Tcp{p) and Tcnt(p), 
the maxima take place at a few po and the fall occurs at the 
densities several times higher. The maximum values of Tc 
range from about 10® K (or even lower) to (2 — 3) x 10^" 
K, depending on the microscopic theoretical model em- 
ployed. The maximum values of Tent are typically lower 
than those of T^p and Tens, due to the weaker nucleon- 
nucleon attraction in the ■^P2 state. 

Instead of studying Tc as a function of p, it is often con- 
venient to consider it versus the nucleon Fermi wavenum- 
ber k = fcHsr = (37r^nN)^^'^, where tin is the number den- 
sity of nucleon species N=n or p. Moreover, instead of 
Tc(fc) one often considers A(fc), the zero-temperature su- 
perfluid gap. For the ^Sq pairing, assuming BCS theory, 
one has A(fc) = Tc(fc)/0.5669. In the case of '^P2 neutron 
pairing the gap depends on the orientation of nucleon mo- 
menta with respect to the quantization axis. Following 
the majority of papers we adopt the pairing with zero 
projection of the total angular momentum on the quanti- 
zation axis. In that case Ant{k) = Tcnt(fc)/0.8416, where 
Ant(fc) is the minimum value of the gap on the neutron 
Fermi surface (correspondin g to t he equator of the Fermi 
sphere; e.g., Yakovlev et al. 1999 ). 

Taking into account a large scatter of Tc{k) provided 
by microscopic theories we do not rely on any particular 
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Fig. 3. Superfluid gaps for models Ip, 2p, and 3p (Fig. 1, 

Table 2) of proton pairing (dot-and-dashed lines) versus 
proton Fermi wavenumber. Solid line 1 and dotted lines 2 
and 3 are the same as in Fig. 2. 
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Fig. 4. Total neutrino emissivity versus density (lower 
horizontal scale) for T = 10^ (solid lines), 3 x 10® (dashes), 
and 10^ K (dot-and-dashes) assuming EOS A at supranu- 
clear densities. Thin lines are for non-superfluid matter; 
thick lines are for matter with Ins, Int, and Ip superflu- 
ids of neutrons and protons. Vertical dotted lines are the 
same as in Fig. 1. Upper horizontal scale shows the radial 
coordinate r in a 1.5 Mq NS. 

microscopic model. Following Papers I and II, we param- 
eterize Tc as 



(fc- 



{k-kof + kf {k-k2f + ki 



(1) 



for ko < k < k2\ and Tc = 0, for k < ko or k > k2- The 
factor To regulates the amplitude of Tc, ko and k2 deter- 
mine positions of the low- and high-density cutoffs, while 
ki and fca specify the shape of Tc(p). All wavenumbers, fc, 
ko, . . . fcs will be expressed in fm~^. We have verified that 
by tuning Tq, fco, •■•^3, our parameterization describes 
accurately numerous results of microscopic calculations. 

In our cooling simulations we consider three models 
of ^Sq proton superfluidity, three models of "'^Sq neutron 
superfluidity, and one model of neutron superfluidity. 
The parameters of the models are given in Table 2, and the 
appropriate Tc(p) are plotted in Fig. ^ We have fco = 
for "'^So pairing. At any given p we choose the neutron 
superfluidity (^So or ■^P2) with higher Tc. 

Models Ins, 2ns, and 3ns of ^Sq neutron superfluid- 
ity are the same as used in Paper II. Models Ins and 2ns 
correspond to about the same, rather strong superfluidity 
(with maximum Tens « 7 x 10^ K). Model 2ns has flat- 
ter maximum and sharper decreasing slopes in the wings 



(near the crust-core interface and the neutron drip point). 
Model 3ns represents a much weaker superfluidity, with 
maximum Tens ~ 2.4 x 10^ K and a narrower density pro- 
file. (One can visualize the radial distributions of Tc in a 
NS by comparing the horizontal scales in Figs. |^ and ^.) 

The superfluid gaps for these models are shown in 
Fig. 1^ versus the neutron Fermi wavenumber. For com- 
parison, we present also three curves provided by micro- 
scopic theories. Solid curve 1 is obtained using BCS the- 
ory with the in-vacuum nn-intcraction (after Lombardo 
& Schulze pOOl| ). This approach yields a very strong su- 
perfluidity with the maximum gap A„s ~ 3 MeV. Dotted 
curves 2 and 3 are calculated using two different mod- 
els of nn-interact ion a ffected by the mediu m polarization 
(Wambach et al. |l993| , Schulze et al. |l996|) . The polariza- 
tion effects strongly reduce the gaps. Microscopic models 
of Ans(fc) are abundant in the literature and the results 
differ considerably (see, e.g.. Fi g. 7 in Lombardo & Schulze 
pool or Fig. 3 in Yakovlev et al. h999|). Our phenomenolog- 
ical curves Ins, 2ns, and 3ns all fall in the range covered by 
theoretical curves. The shapes of the A„si(A;) and Ans3(^) 
curves are typical whereas the decrease of Ans2(fc) with k 
at fc > 1.4 fm~^ is not (too sharp). 

The proton superfluidity curves Ip, 2p, and 3p in Fig. 
1 are similar. The maximum values of Tcp are about 7x 10^ 
K for all three models. Note that model Ip was used in 
Papers I and II. The models differ by the positions of the 
maximum and decreasing slopes of Tcp{p). The decreasing 
slope of model Ip is slightly above the threshold density 
of the direct Urea process (for EOS A), while the slopes 
for models 2p and 3p are shifted to a higher p. The cor- 
responding gaps are shown in Fig. |^. For comparison, in 
Fig. ^ we present the same curves 1, 2, and 3 as in Fig. 
H (the proton gap Aps(fcFp) is expected to be similar to 
the neutron gap Ans(fcFn))- Our models Ip, 2p, and 3p 
are typical for those microscopic theories which adopt a 
moderately strong medium polarization of pp-interaction. 

Finally, the dotted curve in Fig. |l| shows Tcnt{p) for 
our model Int of '^P2 neutron pairing (used in Paper I). 
Microscopic theories give a very large scatter of Tcnt(p), 
and our curve falls within their limits. 

Figure 4 shows the density profile of the neutrino emis- 
sivity at T = 10®, 3x 10®, and 10^ K. Thin fines correspond 
to non-superfluid matter and have three distinct parts. In 
the crust, the emissivity is mainly produced by neutrino 
bremsstrahlung due to the scattering of electrons off nu- 
clei. In the outer core, the emissivity is mainly produced 
by the modified Urea process and is about two orders of 
magnitude higher. In the inner core, it is due to the direct 
Urea process and is higher by another 6-7 magnitudes. 
Thick lines are for superfluid matter (Ins, Int, and Ip su- 
perfluids of neutrons and protons). At T = 10^ K there 
are two large emissivity peaks, near the neutron drip-point 
and the crust-core interface. They are associated with the 
neutrino emission due to ^Sq Cooper pairing of neutrons. 
They are explained by the fact that the Cooper-pairing 
neutrino emissivity is most intense at 0.8 Tc ^ T < Tc, 
and is exponentially small at T ^ Tc (e.g., Yakovlev et al. 
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Table 3. Surface temperatures of eight isolated middle-aged neutron stars inferred from observations 
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Source 



RX J0822-43 

IE 1207-52 

RXJ0002+62 

PSR 0833-45 (Vela) 

PSR 0656+14 

PSR 0633+1748 (Geminga) 

PSR 1055-52 

RX J1856-3754 



[yr] 



3.57 
3.85 
3.95*' 
4.4'=) 

5.00 
5.53 
5.73 
5.95 



[K] 



6.23: 



6 in+°"5 

"•^'-'-0.06 



6.03 



+0.03 
0.03 



+0.02 



5.83 



5 qfi+0.02 

O-yO-0.03 



5.75 

5 

5.72 



0.08 

+0.03 
-0.04 

+0.05 



Model") 



H 
H 
H 
H 

bb 
bb 
bb 



Confid. 
level 



95.5% 
90% 
95.5% 

68% 
90% 

90% 

d) 

d) 



References 



Zavlin et al. (1999) 



Zavlin et al. (1998) 



Zavlin & Pavlov ( |l99S| ) 



Pavlov et al. (2001) 



Possenti et al. (1996) 



Halpern & Wang (1997) 



Ogelman (1995) 



Pons et al. (2001) 



Observations are interpreted either with a hydrogen atmosphere model (H), or with a black body spectrum (bb) 



The mean age taken according to Craig et al. (1997). 



According to Lyne et al. (|l99e|) 



Confidence level is uncertain 
Analytic fit with Si-ash atmosphere model of Pons et al. (2001) 



199E). One can also see a reduction of neutrino emission 
in the outer core by the proton superfluidity. The same 
superfluidity reduces also the direct Urea process near its 
threshold, p ^ pb, but becomes weaker (Fig. ^ and has 
no effect at higher densities. At T = 3 x 10* K the peaks 
associated with ^Sq neutron pairing are weaker, but there 
is a new high peak in the outer core due to ■^P2 neutron 
pairing. At T = 10* K the neutrino emission due to ^Sq 
neutron pairing almost disappears but the emission due to 
^P2 pairing persists. At still lower temperature, the '^P2- 
pairing emissivity in the outer core will have two peaks 
and gradually disappear. Upper horizontal scale gives the 
radial coordinate in a 1.5 Mq NS. The bulk of neutrino 
emission comes evidently from the core ('--^ 11 km in ra- 
dius), and a lower fraction comes from the crust 1 km 
thick). 

The results of cooling calculations are illustrated in 
Figs. and described in Sect. 4. 

3. Observational data 



We will confront theoretical cooling curves with the results 
of observations of thermal emission from eight middle- 
aged isolated NSs. The observational data are the same 
as in Papers I and II. They are summarized in Table 3 
and displayed in Figs. 5, 6, 7, and 10. The three youngest 
objects (RX J0822-43, IE 1207-52, and RX J0002+62) 
are radio-quiet NSs in supernova remnants. The oldest ob- 
ject, RX J1856-3754, is also a radio-quiet NS. The other 
objects. Vela, PSR 0656+14, Geminga, and PSR 1055-52, 
are observed as radio pulsars. The NS ages are either pul- 
sar spindown ages or the estimated supernova ages . The 
age of RX J1856-3754 was estimated by Walter ( ^00 1[ ) 



from the kinematical data (by identifying a possible pre- 
supernova companion in the binary system). We use the 
value t = 9 X 10^ yr mentioned in the subsequent publica- 



tion by Pons et al. (gOOl|). 

For the four youngest sources, the effective surface 
temperatures are obtained from the observed X-ray 
spectra using hydrogen atmosphere models. Such mod- 
els are more consistent with other information on these 
sources (distances, hydrogen column densities, inferred NS 
radii, etc.) than the blackbody model of NS emission. On 
the contrary, for the next three sources we present the val- 
ues of Tg°° inferred using the blackbody spectrum because 
the blackbody model is more consistent for these sources. 
Finally, for RX J1856-3754 we adopt the values inferred 
using the analytic fit with Si-ash atmosphere model of 
Pons et al. (2001). We expect that the large errorbar of 
Tg°° provided by this model reflects poor understanding of 
thermal emission from this source (e.g.. Pons et al. 2001 
Burwitz et al. 2001, Gansicke et al. 2001, Kaplan et al 



pooi| ) 



4. Theory versus observations 

4.1. General remarks 

We have a large scatter of observational limits on for 
the eight sources. Three sources, the youngest RX J0822- 
43, and two oldest, PSR 1055-52 and RX J1856-3754, 
seem to be hot for their ages, while the other ones, espe- 
cially Vela and Geminga, look much colder. Our aim is 
to interpret all observational data with the cooling curves 
using the fixed (the same) EOS and models of the critical 
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temperatures T'c(p) (Sect. 2) for all objects. The results 
are presented in Figs. [5|-|l^. 

In the absence of any superfluidity we would have two 
well-known, distinctly different cooling regimes, slow and 
fast cooling. The slow cooling takes place in low-mass NSs, 
with M < Md. In middle-aged NSs, it goes mainly via 
neutrino emission produced by modified Urea processes. 
For a given EOS in the NS core, the cooling curves of 
middle-aged NSs are almost the same f or all M from about 
1.1 M q to A/d (e.g.. Page & Applegate[l992l Gnedin et al. 



6 6 



2001 ), and they are not very sensitive to EOS. The fast 
cooling occurs if M > M-^ + 0.003 Mq via a very powerful 
direct Urea process (Lattimer et al. 1991). The cooling 



curves are again not too sensitive to the mass and EOS. 
The middle-aged rapidly cooling NSs are much colder than 
the slowly cooling ones. Two examples, for 1.35 and 1.5 
M© non-superfluid NSs (EOS A), are displayed in Fig. 
^ by long dashes. The transition from the slow to fast 
cooling takes place in a very narrow range of M. It is 
demonstrated in Fig. || which shows versus NS mass 
at the age t = 2.5 x ID"* yr of the Vela pulsar. Horizontal 
dotted lines show observational limits on Tg°° for the Vela 
pulsar. One can see a very sharp fall of in the mass 
range Mb < M ^ Mb + 0.003 M© for non-superfluid NSs 
followed by a slow faU at M > Md -t-0.003 Mq. In order to 
explain these observational limits with the non-superfluid 
NS models we should make an unlikely assumption that 
the Vela's mass falls in that narrow mass range. We would 
face the same difficulty with IE 1207-52, RX J0002-h62, 
PSR 0656-1-14, and Geminga. Thus we have five sources 
which exhibit the intermediate case between the slow and 
fast cooling. In the absence of superfluidity, this is highly 
unlikely. 



4.2. Proton superfluidity and the three types of 
cooling neutron stars 

Let us explain the observations (Fig. ||-0) by cooling of 
superfluid NSs. It turns out (Papers I and II) that various 
superfluids affect NS cooling in different ways. Our main 
assumptions would be that (i) the proton superfluidity is 
rather strong a,t p ^ pu, while (ii) the neutron su- 
perfluidity is rather weak (Sect. 4.4). We will discuss the 
superfluid effects step by step starting from the effects of 
proton superfluidity. 

The dot-and-dashed cooling curves in Figs. |^-|^ are 
computed assuming the proton superfluidity alone. We 
adopt the proton pairing Ip in Fig. |^, 2p in Fig. 0, and 3p 
in Fig. 0. We use EOS A in the models in Figs.^ and ^, 
and EOS B in Fig. 



Analyzing Figs. |5|-|§ we see that, generally, the proton 
superfluidity leads to the three cooling regimes (instead 
of two): slow, moderate, and fast. Accordingly, we predict 
three types of cooling NSs with distinctly different proper- 
ties. 




Fig. 5. Observational limits on surface temperatures of 
eight NSs (Table 3) compared with cooling curves for NS 
models (EOS A) with masses from 1.35 to 1.55 M©. Dot- 
and-dashed curves are obtained including proton super- 
fluidity Ip alone. Solid curves include, in addition, model 
Ins of neutron superfluidity. Dotted lines also take into ac- 
count the effect of neutron superfluidity Int. Long-dashed 
lines are for non-superfluid 1.35 and 1.5 M© NSs. 

(I) Low-mass, slowly cooling NSs. The central densi- 
ties pc and masses M — M{pc) of these NSs obey the 
inequalities 



Pc<,pi, M <, Mi^ M{pi); 

the threshold values pi and A/j are specified below. 
(II) Medium-mass, moderately cooling NSs, with 



pi Pc Pih 



Ml <, M <, Mil = M{pii), 



pii and Afn are also specified below. 
(Ill) Massive, rapidly cooling NSs, 



Pi\ ^ Pc 



< 



Mil ^ M < M„ 



(2) 



(3) 



(4) 



where Pcmux and A/max refer to the maximum-mass con- 
figuration (Table 1). 

We will show that the threshold values of pi ji and Mi 
depend on a proton superfluidity model, EOS in the NS 
core, and on a NS age. Let us describe these three cooling 
regimes in more detail. 

(I) We define the slowly cooling NSs as those where 
the direct Urea process is either forbidde n by momentum 
conservation (pc < Pd, Lattimer et al. 1991) or greatly 
suppressed by the strong proton superfluidity. 

In particular, we have the slow cooling for pc < pu and 
M < A/d in the absence of proton superfluidity. This is the 
ordinary slow cooling discussed widely in the literature. It 
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5.5 - 



Fig. 6. Observational limits on surface temperatures of 
NSs compared with cooling curves for NS models (EOS A) 
with several masses M in the presence of proton superflu- 
idity 2p. Dot-and-dashed curves are obtained assuming 
non-superfluid neutrons. Solid curves include, in addition, 
model Ins of neutron superfluidity. neutron pairing is 
neglected. 



6A - 



6.2 




5.6 



Fig. 7. Observational limits on surface temperatures of 
NSs compared with cooling curves for NS models (EOS B) 
with masses from 1.3 to 1.6 M©. Dot-and-dashed curves 
are obtained using model 3p of proton superfluidity. Solid 
curves include, in addition, model Ins of neutron super- 
fluidity. 



is mainly regulated by the neutrino energy losses produced 
by the modified Urea process. 

However, for the conditions displayed in Figs. the 
proton superfluidity is so strong that it almost switches 
off both, the modified Urea process everywhere in the NS 
core and the direct Urea process at p > pu- Then the main 
neutrino emission is produced by neutrino hremsstrahlung 
in neutron-neutron collisions (unaffected by the neutron 
superfluidity in the NS core that is assumed to be weak). 
The hremsstrahlung is less efficient than the modified Urea 
process and leads to an even slower cooling than in a non- 
superfluid NS. We will refer to it as the very slow cooling. 

The analysis shows that, for our proton superfluid 
models, the regime of very slow cooling holds as long as 
the proton critical temperature in the NS center is higher 
than a threshold value: 



rcp(Pc);^Tm(pe). 



(5) 



Comparing the neutrino emissivities of the indicated reac- 



tions (e.g., Yakovlev et al. 1999), we can obtain the simple 
estimates: ripV) ~ 5-5 ^ for p < po, and T^^ {p) - 17T 
for p several per cent higher than pD, where T is the in- 
ternal NS temperature. There is a continuous transition 
between the presented values of Tip'' (p) in the narrow den- 
sity range near pD. 

To make our analysis less abstract we notice that T ^ 
5.5 X 10* K in a very slowly cooling NS at t '-^ 4 x 10^ yr, 
T - 4 X 10* K at t - 2.5 X lO"* yr, and T - 1.5 x 10* K at 
i ~ 4 X 10^ yr. The dependence of Tcp on pc and on the 
NS age is shown in Fig. ||. 

Now we are ready to specify the maximum central den- 
sities Pi and masses Mi of slowly cooling NSs in Eq. (||) . 
For the cases of study, we have pi > pD and Mi > Md 
because the NSs with pc < pd show slow cooling. If 
Pc > Pd, then we may find the density pi on the de- 
creasing, high-density slope of Tcp(p) (Fig. |l|) which corre- 
sponds to Tcp(pi) = Tcp{pi) (Fig. H). It gives the central 
density of the star and Mi ~ M (pi). If pi is formally lower 
than Pd, we set Mi = Md. 

Table 4 shows the values of Mi for EOSs A and B 
and proton superfluids Ip, 2p, and 3p at two NS ages, 
ti = Ax 10^ yr and i2 = 4 x 10^ yr. 

According to Fig. ^, a sufficiently strong proton super- 
fluidity smears out a sharp transition from the slow to fast 
cooling as the mass grows. This effect is illustrated for all 
three proton superfluid models and EOS A at i = 2.5 x 10'' 
yr. In models 2p and 3p the proton superfluidity at the 
direct Urea threshold is very strong. It drastically sup- 
presses the direct Urea process and makes it unimportant. 
In these cases, the direct Urea threshold does not manifest 
a transition to faster cooling. Thus, at t = 2.5 x 10^ yr we 
have Pi pD and Mi ~ Md for proton superfluid Ip, and 
we have pi > pr,, Mi > A/d for superfluids 2p and 3p (Fig. 

I)- 

For the conditions displayed in Figs. the cooling 
curves (dot-and-dashed lines) of all low-mass NSs are very 
similar. For instance, the 1.35 M© curve in Fig. ||is plotted 
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Table 4. Masses Miji which separate slow, moderate, and 
rapid cooling models for ti = 4 x 10^ and t2 = 4 x 10^ yr 



EOS 


Proton 


Mi/Mq 


Mii/Mq 




pairing 


tl 


t2 


tl t2 


A 


Ip 


Mn 


1.4 


1.52 1.53 


A 


2p 


Mo 


1.55 


1.64 1.64 


A 


3p 


Mo 


1.77 


1.83 1.84 


B 


Ip 


Mo 


Mo 


Mo Mo 


B 


2p 


Mo 


Mo 


Mo Mo 


B 


3p 


Mo 


1.55 


1.62 1.62 



just as an example; all the curves are almost identical in 
the mass range 1.1 Mq ^ M < Mj. Moreover, the curves 
are not too sensitive to EOS and are insensitive to the 
exact values of the proton critical temperature Tcp, as long 
as the inequality (||) holds. They are noticeably higher 
than the analogous cooling curves of the ordinary slow 
cooling in the absence of superfluidity (e.g.. Fig. ^). 

For the conditions in Figs. (as in Papers I and 
II) we may explain the three relatively hot sources, RX 
J0822-43, PSR 1055-52, and RX J1856-3754, by these 
very-slow-cooling models with a strong proton superfluid- 
ity. Thus, we assume that the indicated sources are low- 
mass NSs. We discuss this explanation further in Sect. 
4.3. 

(II) We define the moderately cooling stars as the NSs 
which possess central kernels where the direct Urea pro- 
cess is allowed but moderately suppressed by proton su- 
perfluidity. The existence of a representative class of these 
NSs is solely due to proton superfluidity. 

Our analysis shows that, for our cooling models, the 
proton critical temperature in the center of a medium- 
mass NS should roughly satisfy the inequality 



Tm(pe);$Tep(pe);$Tm(pe), 



(6) 



with ri"' - ST. Thus, we may introduce pn which cor- 
responds to rcp(pii) = Tc^\pii) and determines Afn, the 
maximum mass of moderately cooling NSs in Eq. (^ . The 
values of depend mainly on a given EOS and slightly 
on a NS age (Fig. ||). The ranges of mass and density in 
Eq. (^ depend also on a model of Tcp(p). 

The values of M\\ are also given in Table 4, along with 
Ml. For EOS B, the critical temperature Tcp{p) of Ip and 
2p proton superfluids vanishes at p ^ pD • Then we have a 
sharp transition from the slow to fast cooling in a narrow 
mass range just as in the absence of the superfluidity (Sect. 
4.1; Ml « Mil ~ Mu), and the regime of moderate cooling 
is almost absent. 

The surface temperatures of the medium-mass (mod- 
erately cooling) NSs are governed by proton superfluidity 
in the NS central kernels, p ^ pi. One can observe (Figs. 
^-^) a steady decrease of surface temperatures with in- 
creasing M. If we fix the proton superfluidity and EOS 
(provided they allow for the moderate cooling) we can de- 
termine (Papers I and II) the mass of any moderately cool- 
ing NS, which means "weighing" NSs. In this fashion we 



can weigh five isolated NSs (IE 1207-52, RX J0002-h62, 
Vela, PSR 0656-f 14, and Gcminga) using either EOS A 
and the proton superfluids Ip, 2p, and 3p, or EOS B 
and the superfluid 3p. Thus, we assume that the indicated 
sources are moderately cooling NSs. For instance, adopt- 
ing EOS A and proton superfluid Ip (Fig. |^) we obtain 
the masses in the range from « 1.36 M© (for IE 1207-52) 
to « 1.465 Mq (for Vela and Geminga). For EOS A and 
proton superfluid 2p (Fig. ^) we naturally obtain higher 
masses of the same sources. Obviously, the properties of 
moderately cooling NSs are extremely sensitive to the de- 
creasing slope of Tcp(p) in the temperature range from 



Tcp^ to Tcp (Fig. H), or in the density range from pi to 
pii (and insensitive to the details of Tcp(p) outside this 
range) . 

(Ill) Massive NSs show fast cooling similar to the fast 
cooling of non-superfluid NSs. These stars have central 
kernels where the direct Urea process is either unaffected 
or weakly suppressed by the proton superfluidity. In such 
kernels, Tcp(p) ^ Tcl^\ The central densities and masses 
of rapidly cooling NSs lie in the range given by Eq. (^. 
The thermal evolution of rapidly cooling NSs is not very 
sensitive to the model of T'cp(p) and to EOS in the stellar 
core. Note that if pcmax < Pii, the rapidly coohng NSs 
do not exist. In the frame of our interpretation, no NS 
observed so far can be assigned to this class. 



4.3. Crustal superfluidity and slow cooling 

As the next step, we retain proton superfluidity and add 
^So neutron superfluidity Ins in the NS crust and outer- 
most core. In this case we obtain the solid curves in Figs. 
Ij-Q. For the moderately or rapidly cooling middle-aged 
NSs {M > Ml) they are fairly close to the dot-and-dashed 
curves. This is quite expected (e.g., Gnedin et al. 2001): 
the ^So neutron superfluidity is mainly located in the NS 
crust which is much less massive than the NS core. Thus, 
the crustal superfluidity does not affect noticeably our 
interpretation of IE 1207-43, RX J0002-f62, Vela, PSR 
0656-1-14, and Geminga in terms of moderately cooling 
NSs. 

However, as pointed out in Paper II, this crustal super- 
fluidity strongly affects the slow cooling of low-mass NSs 
(M < Ml). Its effects are twofold. First, at t ^3x10^ yr 
the neutrino luminosity due to ^Sq pairing of neutrons may 
dominate the sufficiently low neutrino luminosity of the 
stellar core (see Fig. ^ and Paper II, for details). Second, at 
t ^ 10^ yr the ^Sg neutron superfluidity reduces the heat 
capacity of the crust. Both effects accelerate NS cooling 
and decrease (Figs. ||-0) violating our interpretation 
of the three suflfciently hot sources, RX J0822-43, PSR 
1055-52, and RX J1856-3754. The interpretation of RX 
J1856-3754 is affected to a lesser extent, as a consequence 
of the rather large errorbar of for this source (Sect. 
3). 

Let us demonstrate that our interpretation can be res- 
cued by the appropriate choice of Tcnsif)- For this pur- 
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pose we focus on the interpretation of RX J0822-43, PSR 
1055-52, and RX J1856-3754, as the very slowly cooling 
NSs (M < Ml). 

For certainty, let us take EOS B, M = 1.3 Mq, and 
proton superfluid 3p. The results are presented in Fig. ^ 
(cf with Fig. 3 of Paper II). The dot-and-dashed line is 
the same as in Fig. ^ and neglects the crustal neutron su- 
perfluidity. Thick solid line is also the same as in Fig. |7| It 
includes an additional effect of crustal superfluid Ins and 
lies below the observational limits on for the sources 
in question (or almost below in case of RX J1856-3754). 
To keep the proposed interpretation of the three sources 
we must raise the cooling curves calculated including the 
crustal superfluidity. To this aim, we must suppress the 
neutrino emission associated with ^Sq pairing of neutrons 
(Fig. U). Recall that in a middle-aged NS this emission is 
mainly generated (Fig. ^ in two relatively narrow layers, 
near the neutron drip point and near the crust-boundary 
interface, where the local NS temperature T is just below 
Tcnsip)- Since the Cooper-pairing neutrino luminosity is 
roughly proportional to the widths of these emitting lay- 
ers, we can achieve our goal by reducing their widths. This 
can be done by setting T^"^^^ higher and by making Tcns{p) 
decrease sharper in the wings (see Paper II, for details). 

For example, taking crustal superfluid 2ns instead of 
Ins (Figs. |l| and ||) we obtain the dashed cooling curve 
in Fig. |l^ which comes much closer to the dot-and-dashed 
curve than the thick solid curve (model Ins). (Another ex- 
ample: shifting Tcns(p) for model 2ns into the crust would 
additionally raise the curve towards the dot-and-dashed 
one.) Note that the cooling curves are insensitive to the 
details of Tcns(p) profile near the maximum, as long as 
^cns'' ^ 5 X 10^ K, but they are extremely sensitive to the 
decreasing slopes of Tcns{p)- On the other hand, by taking 
the smoother and lower Tcnsip), model 3ns, we obtain a 
colder NS than needed for the interpretation of observa- 
tions (long-dash line in Fig. |l^. Therefore, ^Sq neutron 
superfluidity with maximum Tens'' < 5 x 10^ K and/or 
with smoothly decreasing slopes of the Tens (p) proflle near 
the crust-core interface and neutron drip point violates the 
proposed interpretation of the observational data. 

Let us stress that the observations of RX J0822-43, 
PSR 1055-52, and RX J1856-3754 can be fitted even with 
our initial model Ins of the crustal neutron superfluidity. 
The high surface temperature of RX J0822-43 can be ex- 
plained assuming additionally the presence of a low-mass 
(2 X 10~^^ Mq) heat-blanketing surface envelope of hydro- 
gen or helium. This effect is modeled using the results of 
Potekhin et al. (1997). Light elements increase the elec- 
tron thermal conductivity of NS surface layers and raise 
at the neutrino cooling stage (curve acc in Fig. p^ . 
In order to explain the observations of PSR 1055-52 and 
RX J1856~3754, we can assume again model Ins of crustal 
superfluidity, iron surface and the dipole surface magnetic 
field (~ 10^^ G at the magnetic pole; line mag in Fig. 
[l0| ). Such a field makes the NS surface layers overall less 
heat-transparent (Potekhin & Yakovlev 2001), rising 
at i ^ 3 X 10^ yr. Note that the dipole field ^ lO^^ q 



has the opposite effect, resembling the effect of the sur- 
face envelope of light elements. Therefore, we can addi- 
tionally vary cooling curves by assuming the presence of 
light elements and/or the magnetic field on the NS sur- 
face. However, these variations are less pronounced than 
those due to nucleon superfiuidity. For instance, we cannot 
reconcile the cooling curves with the observations of PSR 
1055-52 assuming model 3ns of the crustal superfluidity 
with any surface magnetic field. 

4.4. pairing of neutrons in tlie NS core 

Now we focus on the effects of '^P2 neutron pairing, which 
were neglected so far. They are illustrated in Fig. ||, as 
an example. They would be qualitatively similar for the 
other cooling models in Figs. ^ and ^ In Fig. |^ we take the 
cooling models obtained including proton superfluidity Ip 
and crustal superfluidity Ins and add the ''P2 neutron su- 
perfluidity (model Int, Table 2) in the core. We have the 
same (solid) cooling curves for the young NSs which have 
the internal temperatures T above the maximum value of 
T™f^ « 3 X 10*^ K. However, when T falls below T^f'^, we 
obtain (dots) a strong acceleration of the cooling associ- 
ated with the powerful neutrino emission due to ^P2 neu- 
tron pairing (Fig. ji]). This emission complicates our inter- 
pretation of older sources, PSR 0656-1-14, Geminga, PSR 
1055-52, and RX J1856-3754. It may induce really fast 
cooling of such sources even if their mass is low, M < Md 
(Sect. 4.2). To avoid this difficulty we assume (Paper I) 
weak '^P2 pairing, Tcnt{p), with maximum T^f^ < 10^ K. 
Then, it does not affect the proposed interpretation. 

5. Summary 

Let us summarize the effects of the three types of super- 
fluids on NS cooling: 

(a) Strong proton superfluidity in the NS cores, com- 
bined with the direct Urea process at p > pu, separates 
the cooling models into three types: (I) slowly cooling, low- 
mass NSs (M ^ -Mi); (II) moderately cooling, medium- 
mass NSs (Mi < M <, Mil); (HI) rapidly coohng, mas- 
sive NSs (M Mil). These models have distinctly differ- 
ent properties (Sect. 4.2). The regime of moderate cooling 
cannot be realized without the proton superfluidity. 

(b) Strong proton superfluidity in the NS core is re- 
quired to interpret the observational data on the three 
sources, RX J0822-43, PSR 1055-52, and RX J1856-3754, 
hot for their ages, as the very slowly coohng NSs (Sects. 
4.2, 4.3). Within this interpretation, all three sources may 
have masses from about 1.1 M© to Mi; it would be diffi- 
cult to determine their masses exactly or distinguish EOS 
in the NS core from the cooling models. 

(c) Strong proton superfluidity is needed to inter- 
pret observations of the other sources, IE 1207-52, RX 
J0002-h62, Vela, PSR 0656-hl4, and Geminga, as the 
medium-mass NSs. This allows one to "weigh" these NSs, 
i.e., determine their masses, for a given model of Tcp{p). 
The weighing is very sensitive to the decreasing slope of 
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Fig. 8. Surface temperatures of NS models (EOS A) at 
t = 2.5 X 10'* yr, the age of the Vela pulsar, versus stel- 
lar mass M. The curves are calculated (from bottom to 
top) assuming either no superfluidity, or proton superflu- 
idities Ip, 2p, 3p (and neglecting neutron superfluidity). 
The vertical dotted line shows the threshold mass of 
opening the direct Urea process. Horizontal dotted lines 
are observational limits on the Vela's surface temperature 
(Table 3). 



Tcp{p) in the density range pi ^ p pii, and it depends 
also on the EOS in the NS core (Sect. 4.2). 

(d) Strong or moderate '^P2 neutron superfluidity in 
the NS core initiates rapid cooling due to the neutrino 
emission resulted from neutron pairing. This invalidates 
the proposed interpretation of the old sources like PSR 
0656-hl4, Geminga, PSR 1055-52, and RX J1856-3754. 
To save the interpretation, we assume a weak ■^P2 neutron 
superfluidity, T^f < 10^ K (Sect. 4.4). 

(e) *So neutron superfluidity in the crust may initi- 
ate a strong Cooper-pairing neutrino emission, decrease 
substantially of the slowly cooling NSs, and weaken 
our interpretation of RX J0822-43, PSR 1055-52, and RX 
J1856-3754 (although it does not affect significantly the 
moderate or fast cooling). We can save the interpretation 
by assuming that the maximum of the critical tempera- 
ture profile Tcnsip) is not too small (T™^'^ ^ 5 x 10^ K) 
and the profile decreases sharply in the wings (Sect. 4.3). 

(f) The interpretation of the slowly cooling sources is 
sensitive to the presence of the surface magnetic fields 
and/or heat-blanketing surface layer composed of light el- 
ements (Sect. 4.3). 

(g) No isolated middle-aged NSs observed so far can 
be identified as a rapidly cooling NS. In the frame of our 
models, these NSs do not exist for those EOSs and super- 
fluid T'cp(p) for which M^ax < -^^ii- 



Fig. 9. Two characteristic critical temperatures of proton 
superfluidity, Tcp and Tcp \ versus NS central density pc 
for EOS A and three NS ages, 4 x 10^, 2.5 x 10", and 
4 x 10^ yr. Also shown are Tcp[p) (dot-and-dashed) curves 
for proton superfluids Ip, 2p, and 3p. Their intersections 
with Tcp'^^^ at Pc > Pd separates the low-, medium-, and 
high-mass NS models (see text for details). 



If our interpretation is correct, we can make the fol- 
lowing conclusions on the properties of dense matter in 
NS interiors. 

(i) Strong proton superfluidity we need is in favor of a 
not too large symmetry energy at supranuclear densities 
(Paper I). A very large symmetry energy would mean a 
high proton fraction which would suppress proton pairing. 
On the other hand, the symmetry energy should not be 
too small to open the direct Urea process at p > pD- 

(ii) Weak neutron pairing is in favor of a not too 
soft EOS in the NS core (Paper I). The softness would 
mean a strong attractive nn interaction and, therefore, 
strong neutron pairing. 

(iii) Specific features of the crustal neutron superflu- 
idity we adopt are in favor of those microscopic theories 
which predict Tc„s(p) profiles with T^^f" ^ 5 x 10^ K (or 
Ac„(A:) profiles with A^f'^ ^ 1 MeV, see Fig. g). This is in 
line with many microscopic calculations of the superfluid 
gaps which include the medium polarization effects in nn 
interaction (e.g., Lombardo & Schulze 2001). However, the 
reduction of the gap by the medium polarization should 
not be too strong, and the decreasing slope of A(fc) should 
be rather sharp. These requirements constrain the micro- 
scopic theories. 

The proposed interpretation of the observations relates 
the inferred NS masses to the superfluid properties of NS 
interiors. By varying EOS and the proton critical tempera- 
ture, we can attribute different masses to the same sources. 
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Fig. 10. Cooling curves of 1.3 NS model versus ob- 
servations of RX J0822-43, PSR 1055-52, and RX J1856" 
3754. Dot-and-dashed curve: proton superfluidity 3p in 
the NS core. Solid, short-dashed, and long-dashed curves 
include, in addition, models Ins, 2ns, and 3ns of crustal 
neutron superfluidity, respectively. Thick solid line is the 
same as m Fig. 0. Thin sohd curve acc is calculated as- 
suming the presence of 2 x 10""'^"'^ M0 of hydrogen on the 
NS surface. Thin solid curve mag is obtained assuming 
the dipole surface magnetic field (10^^ G at the magnetic 
pole). 



If, on the other hand, we knew the range of masses of the 
cooling middle-aged NSs we would be able to draw defi- 
nite conclusions on the superfluid state of their interiors, 
flrst of all, on the proton critical temperature, Tcp(p). 

Our analysis may seem too simplified because we ne- 
glect a possible presence of other particles in the NS cores 
(muons, hyperons, quarks). We expect that the inclusion 
of other particles and the effects of superfluidity of hyper- 
ons or quarks will complicate theoretical analysis but will 
not change our basic conclusion on the existence of the 
slowly, moderately, and rapidly cooling NSs. 

Our calculations show that the cooling of middle-aged 
NSs with M < Ml is sensitive to the density profile of free 
neutrons near the crust bottom and neutron drip point. 
We have used only one model of the free-neutron distribu- 
tion in the crust, assuming the atomic nuclei to be spheri- 
cal at the crust bottom. It would be interesting to consider 
the models of crust matter with non-spherical nuclei (e.g., 
Pethick & Ravenhall 1995 ) and the effects of superfluidity 
of nucleons conflned in the atomic nuclei in the NS crust. 

Let us stress that determination of from observa- 
tional data is a very complicated problem (as described in 
part by Yakovlev et al. 1999 ). It requires very high-quality 
data and theoretical models of NS atmospheres. Thus, the 



current values of may change substantially after the 
forthcoming observations and new theoretical modeling. 
These changes may affect our interpretation of the obser- 
vational data, first of aU, of RX J0822-43, PSR 1055-52, 
and RX J1856-3754. For instance, RX J1856-3754 may 
have a colder surface {T^ ~ 0.25 MK), than assumed in 



the above analysis, with a hot spot (e.g.. Pons et al. 2001 
Burwitz et al. pOOl , Gansicke et al. 2001). If confirmed 
the lower might be explained by the effect of ■^P2 neu- 
tron pairing (Fig. ^). We expect that future observations 
of the thermal emission from these sources will be crucial 
for understanding the superfiuid properties of NS matter. 
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